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Stochastic  Differential  Equations  for  Neuronal  Behaviour 


In  this  article  we  extend  the  recent  work  of  Kallianpur  and  Wolpert 
[5]  modeling  the  behavior  of  neurons  by  means  of  stochastic 
differential  equations  on  the  dual  of  a  nuclear  space.  The 
extensions  will  cover  nuclear  spaces  of  a  more  general  structure  and 
will  apply  to  models  described  in  terms  of  more  general  differential 
operators.  A  second  objective  of  our  work  is  to  present  a 
theoretical  framework  which  will  include  the  model  recently  proposed 
and  heuristically  investigated  by  Wan  &  Tuckwell  in  [9].  We 
illustrate  our  approach  and  its  application  by  giving  a  rigorous 
treatment  of  the  Wan  &  Tuckwell  model.  But  first  we  shall  briefly 
describe  the  neurophysiological  context.  For  a  more  detailed  account, 
we  refer  to  [5]  and  the  references  therein.  In  our  description  we 
shAll  follow  the  introduction  in  [5]. 


A  neuron  is  a  cell  whose  principal  function  is  to  transmit 
information  along  its  considerable  length,  which  often  exceeds  one 
meter.  "Information"  is  represented  by  changing  amplitudes  of 
electrical  voltage  potentials  across  the  cell  wall.  A  quiescent 
neuron  will  exhibit  a  resting  potential  of  about  60  mV,  the  inside 
more  negative  than  the  outside.  Under  certain  circumstances  the 
voltage  potential  in  the  neuron  dendrite  will  rise  above  a  threshold 
point  at  which  positive  feedback  causes  a  pulse  of  up  to  100  mV  to 
appear  at  the  base  of  the  dendrite;  this  pulse  is  transmitted  rapidly 


along  the  body  and  down  the  axon  of  the  cell  until  it  reaches  the 
so-called  "pre— synaptic  terminals"  at  the  other  end  of  the  neuron. 
Here  the  pulse  causes  tiny  vesicles  filled  with  chemicals  called 
"neurotransmitters”  to  empty  into  the  narrow  gaps  between  the 
presynaptic  terminals  and  the  dentrites  of  other  neurons.  When  these 
chemicals  diffuse  across  the  gap  and  hit  the  neighboring  neurons' 
dendrites,  they  may  cause  the  voltage  potential  in  these  dentrites  to 
rise  above  a  threshold  point  and  initiate  another  pulse. 


Let  ^»(t,x)  represent  the  difference  between  the  voltage  potential  at 
time  t  at  the  location  x  e  X  (■  surface  of  the  neuron)  and  the 
resting  potential  of  about  —60  mV.  As  time  passes,  5  evolves  due  to 
two  separate  causes: 

(i)  Diffusion  and  leaks:  Depending  on  the  nature  of  X,  the 

electrical  properties  of  the  cell  wall  may  be  approximated  by 

postulating  a  contraction  semigroup  {Tfc>  on  L^(X,7*)  where  ?  is  a 

/ 

suitable  6— finite  measure  on  X.  For  example,  if  X  ■  [0,b],  core 
conductor  theory  suggest  the  semigroup  corresponding  to  the  diffusion 
equation 


-*St  +  <8,4  >  0) 


with  Neumann  (or  insulating)  boundary  conditions  at  both  ends.  In 
neural  material  like  heart  muscle  in  which  electrical  signals  can 
travel  more  easily  in  some  directions  than  in  others,  the  Laplacian 
should  be  replaced  by  a  more  general  second— order  elliptic  operator. 


(ii)  Random  fluctuations:  Every  now  and  then  a  burst  of 


neurotransmitter  will  hit  some  place  or  another  on  the  membrane  and 
suddenly  the  membrane  potential  will  jump  up  or  down  by  a  random 
amount  at  a  random  time  and  location.  It  is  believed  that  these 
random  jumps  are  quite  small  and  quite  frequent,  making  it  reasonable 
to  hope  that  they  can  be  modeled  by  a  Gaussian  noise  process;  in  any 
case  the  arrivals  at  distant  locations  or  in  disjoint  time  intervals 
are  believed  to  be  aproximately  independent,  justifying  their 
modeling  as  a  mixture  of  Poisson  processes  or  as  a  generalised 
Poisson  process. 

Taking  into  account  these  random  fluctuations  one  arrives  at  the 
following  stochastic  partial  differential  equation  (SPDB)  for  : 

d£(t,x)  «  A^(t,x)dt  +  dMt;  initial  condition 

Here  A  is  a  suitable  partial  differential  operator  in  spatial 
I  2 

coordinates  and  Mt  is  an  L  -semimartingale.  However,  even  for  very 
simple  choices  of  A  (e.g.  A  ■  I_A  *n  two  dimensions;  see  [8])  a 
solution  may  exist  only  in  the  form  of  a  generalized  stochastic 
process  i.e.  a  process  taking  values  in  the  dual  of  a  space  $  of 

"testf unctions".  The  relevant  space  of  testfunctions  can  usually  not 

\ 

be  assumed  to  be  the  Schwartz  space  of  all  infinitely  differentiable 
rapidly  decreasing  functions  (see  e.g.  [5])  and  therefore  we  shall 
take  $  to  be  a  general  countably  Hilbert  nuclear  space.  The  linear 
SPDEs  appropriate  for  this  purpose  have  been  investigated  in  [2] 
where  an  existence  and  uniqueness  result  is  given  for  equations 
driven  by  a  semimartingale  on  . 


4 


In  [5]  a  restricted  class  of  differential  operators  was  considered; 
namely  those  which  generate  a  self adjoint  contraction  semigroup  whose 
resolvent  has  a  power  which  is  Hilbert-Schmidt.  In  this  case  there  is 
a  canonical  nuclear  space  upon  which  the  SPDE  has  a  very  manageable 
form.  However ,  the  structure  of  the  nuclear  space  is  completely 
de^rmined  by  the  operator  A,  and  therefore  we  shall  present  general 
results  which  are  independent  of  the  structure  of  $  and  which  will 
also  permit  a  wider  class  of  differential  operators  to  be  considered. 

Kallianpur  and  Wolpert  ([5])  used  a  Poisson  process  N(AxBx(0,t])  to 
represent  the  number  of  voltage  pulses  of  size  a  6  A  arriving  at 
sites  x  e  B  C.  X  («  surface  of  the  neuron)  at  times  prior  to  t. 

Here,  we  adopt  the  point  of  view  that,  in  practice,  one  can  only 
"average"  over  the  sites.  Therefore  it  seems  more  realistic  to  assume 

that, the  arrival  sites  are  given  by  "generalized  functions" 

/ 

(distributions)  e  A  ^  »  rather  than  by  points  x  on  the  surface 

of  the  neuron  membrane  X.  As  we  shall  see,  this  approach  will  also 
offer  the  advantage  of  enlarging  the  class  of  possible  models. 

To  pursue  this  idea  let  us  consider  a  real  rigged  Hilbert  space 
(see  Gel'fand  &  Vilenkin  [3]  p  for  definition).  Let 
fl(<£' )  denote  the  Borel  field  on  and  recall  that  B(<jjr  >  is  the 
same  whether  we  use  the  weakly  or  the  strongly  open  sets  in  to 
define  it. 

-To  avoid  possible  confusion  with  inner  products  we  shall.  adopt  the 
notation  that  for  ^  €  $  and  rp  6  <$' ,  will  denote  the  value  of 


the  functional  evaluated  at 


rLet  A  6  8($'>  and*  let,  for  each  n  6  ■,  nn  be  a  C-f inite  positive 
measure  on  (RxA*S(R)x8(A) )  satisfying: 

The  mapping:  Qn  :  fixjj  -»  R  defined  by 

Qn<4r1»)  *  a2J}[$]  l^n|r]yn(dadrj)  is  continuous  on 

Let  Nn  be  a  Poisson  random  measure  on 

(RxAx[0,oo);8(R)x8(A)xB([0,oo) )  with  intensity  measure  |in(dad]^)dt 
(a  6  R,  6  A#  t  e  [0,oo))  (see  e.g.  Ikeda  and  Watanabe  [4],  page 
42). 

Let  Nn(dad]^ds)  »  Nn(dadJ^ds)  —  nn(dadl}  )ds  and  put 

Let  m11  e  ,  and  define  x£(<f>)  *  tmn[<j>]  +  y£(4>),-  $  e  $• 

Then,  for  each  ^  6  $>,  x"(<J>)  is  a  real  CADLAG  semimartingale  satisfying 
E<x” (4») ) 2  -  t2mn(4>]2  +  tQn(«j>,<^). 

Since  Qn  is  continuous  on  $x$,  the  Kernel  theorem  for  nuclear  spaces 
(see  Gel' f and  &  Vilenkin  [3],  page  74)  yields  the  existence  of 
r(n)  e  ■  and  C(n)  >  0  such  that 


mn[4>]2  +  Qn«M)  <  C(n)|J«H|2(n), 


V  <f>  e 


We  shall  henceforth  assume  that  the  same  r  and  C  will  do  for  all 
n  6  ■,  i.e.  we  suppose  that  there  exists  r2  e  ■  ,  C  >  0  such  that 

(1)  mn[4>J2  +  Qn(<|>,<fr)  <  c||M2  V  ne  i  y  4>  e  f . 

2 

From  Doob's  Martingale  inequality  we  deduce  that,  for  any  T  >  0, 

Esup  <x”<4>))2  <  2C(  4T  +  2t2)||^||2  yne*  v<j»ei 
o<t<T  r2 

and  therefore  Theorem  III. 1.12  and  Remark  7  of  [2]  yields  the 

existence  of  q  e  I,  q  >  r^  (independent  of  n)  and  a  $  valued 

CADLAG  L2 -semimartingale  x"  satisfying  x"t$]  ■  x£($)  V  t>0  (a.s.)  V 

♦  «  I-  X?,T  >-  <x?>teto,Ti'  T  >  «• 

Let  m  e  and  let  Q  :  Jx^  4  I  be  a  continuous  bilinear  symmetric 

/ 

functional  satisfying 

(2)  m[$]2  +  Q ( 4> , 4> )  <  C||<H|2  . 

r2 

It  can  then  be  shown  that  there  exists  a  ^-valued  process  with 

\ 

independent  increments  and  characteristic  functional  given  by 
exp(  itm[4>]  -  t/2Q($,$)).  This  was  shown  by  V.  Perez-Abreu  [7]  for 
the  case  m«0  and  a  nuclear  space  of  a  special  structure.  The  general 
case  may  be  deduced  from  theorem  III. 1.12  in  [2].  We  shall  henceforth 
call  w  a  -valued  Wiener  process  with  parameters  m  and  Q. 

It  can  be  shown  from  (1)  and  (2)  (See  [2]  theorem  III. 1.12  and  remark 
7 )  that  we  may  choose  q  >  r,  such  that 


Xn'T  6  pap,?],*!)  _)  P-a.s.  V  n  e  ■  V  T  >  0  and 

4 

WT  e  C([0,T],3  >  P-a.s.  v  T  >  0. 

4 

Let  p”  denote  the  measure  induced  on  D(  [ 0,T]  ,<|_g)  by  Xn'T  and  let  PT 
denote  the  measure  induced  on  C( [0,T]  »$_q)  ^  D([0,T]  ,<£_q)  by  WT . 


Let  Qnfmn  satisfy  (1).  Then,  for  every  T>  0,  the  family  {P®  :  n  e  *> 
is  tight  on  D<  [Q,T]  ,<ft  ) . 

PROOF: 


By  Mitoma,  [6],  theorem  5.3.2.  and  remark  Rl,  (page  996/7) it  is 
sufficient  to  show  that 

i 

/ 

(a)  V  <M  I  :  {Xn'T[$]  s  n  e  ■)  is  tight  on  D([0,T],R) 


and 

(b)  3k>q:Ve>0Vp>03^>0! 

Hit  <  4  »>  p(sup  |xj[^]|  >e><p  V  n  e  ■ 

°<t<T  1 

For  part  (a),  by  Billingsley  11]  theorem  15.3  page  125,  it  is 
sufficient  to  show  that  V  $  e  $! 


(ai ) 


Vn>03a>0:  P  ( sup  |  xjH>]  |  >  a)  <  p 

1  0<t<T  C  1 

r  (aii )  V  6  >  0,  J|  >  0  3  ^  e  (0,T)  3  nQ  e  ■: 

P { sup  min{|xjl4>]  -  xj  l$]|,|x£  -  xj[<|>]|}  >  e>  <  p 

ti<t<t2  *  zl  z2  1 

and 

P ( sup  ,  |x"[4>l  -  x"[<i>l|  >  e)  <  r)  V  n  >  n 

s,te[o,6)  St  (  O 

and 


P(sup  .  |x"[4»]  “  xj[4>]|  >  e)  <  D  V  n  >  n0 

s,te[T-d>,T)  sc  < 

Fix  <i>  e  $r  and  let  >  0,  e  >  0.  Then, 

)  1  n  2 

P( sup  |xj[«j)]|  >  a)  <  -5  E(sup  |X"[<|>]|)Z 
teto,Ti  c  a*  te[o,T] 

2 

<  -r  E( sup  ( t2mnt 4»3 2  +  Yj[^>]2)) 

a *  te[0#T] 

2C 

<  -r(T2mn[«|>]2  +  4TQn (<f>,<f>)) 
a^ 

2  n  o 

<  -x(Tz  +  4T)  C||<|>||;  V  n  e  ■  (by  1) 

a4  r2 

2 

<  D  for  a2  >  — <T2  +  4T)C|j$||  2 

*  I?  r2 

Next,  let  D:-{t>0:  t^  <  t  <  t^  and  tj-t^  -  Then 


V  n€H 


V  n  >nG 


P{sup  min{|X?[4>l  —  t4>l|»|X?  [^1  —  x”t^l|)  >.  6  } 

teo  t  fcl  2 

<<  e"2  e{  sup  min{|xj[4>l  -  x?  t4>l  I » IxJ  [+]  -  x£[$]  |  }>2 
teo  12 

<  e“2  E{  sup  min{|xj?[4>]  —  X?  I^]|2f|x”  [$]  —  x”[^]|2) 

tSD  t  rl  '  2 

<  e"2  min (Esup  |x"  C$]|2,Esup  |X?_.  [ ^3 1 2 ) 

teD  teD  a. 


2 

<  — 7inin{sup  ((t  —  t,)2mn[$]2)  +  Esup  |Yj?_t  1 4* 3  | ^ r 
e2  t€D  1  t€D  1 

sup  ( (t,  —  t)2mn[$]2)  +  Esup  |Y?  _ 1 1 4>  1  1^3 

teD  2  teD  z2 


=  A(^2mn[4>]2  +  4^Qn(<t>,<|))) 

6  i 
2 

<  (&2  +  4i)c||  <j>  ||  2  V  n  >  1  (by  <1>> 

"  £2  r2 

2 

<  n  v  n  e  ■  if  (i2  +  4^>)  <  <— ~ — c|M>||2  j-"1. 

(  e^  r2 


Further , 


P ( sup  |x"[«t>l  -  x”t4>3 1  >  6) 

sfteio,6) 

<  -^E(sup  .  | x”l 4>3  -  x"c 4) 3  > 2 

“  sr  s,te[0,a) 

X  2 

<  — x2  (£2mn[  $3 2  +  4^Qn(4>,<|)) )  <  5 ( +  4<^)C |(  4»  ||  ^  V  n  >  1 

-  or  s  2 


10 


2 

<  n  V  n  e  ■  if  i2  +  4^  <  ( — ~CH\\1  )“x. 

rje^  r2 


Similarly, 


P(sup  |x"[<|>]  -  X^I4>]  J  >  6) 

s,te[T-d,Tl  8  t  r  - 

<  2e“2(i2mn[^]2  +  4^Qn<4>,4>) )  <  2e”2(i2  +  4<^)C ||  4> ||  2  V  n  >  1 

r <%  ~ * 


<  1}  V  n  e  ■  if  i2  +  4<&  <  (—5 — C II  4» II  ?  J*1. 

6*  r2 


Hence  (ai)  and  (aii)  are  satisfied  for 

A2  +  <A  <  (2IJ-1e_2cH#2  >-1,  and  n  -  1.  This 


proves  (a). 


<b) 


s  Fix  4  e  $  and  let  e,  J}  >  0. 


Then  P< sup  |xj[4>]|  >  e)  <  e”2E(sup  |x?[d>]|): 
tetofT]  te[o,T]  t 

<  e~22(T2mn[<|>J2  +  4TQn<4»,<|>) ) 


<  e  22<T2  +  4T)C||<|>||22  <  e  22CT2  +  4T)c||<|>||2  (by  <i)> 


<  1]  V  n  >  1  if 


q  -  A2 


2 (T  +4T)C 


This  completes  the  proof  of  Proposition  1.0. 


I 


1.1.  THEOREM! 


11 


Suppose  that,  in  addition  to  assumption  (1), 


Qn (4>,<)»)  — ¥  q<4>,$)  v  <►  e 

n-»co 

lim^  /^lad^I^J)  |3^n(dad]^)  -  0  V  i  e 


mnt4>] - >  mt$]  V  ♦  6|. 

n-»oo 


Then,  for  any  T  >  0,  we  have 


p“  =s\  p 

n-»oo 

PROOF: 


Fix  T  >  0.  Let  tj  <  tj  < 
K  6  «  fixed. 


We  must  show  that 


<  tR  e  [ 0,T]  and  ^1,...4k  e  $  for 


(ii) 


(x”  converges  in  distribution  to  (Wfc  1^1) 

^  Jc 


{P^  :  n  e  ■>  is  tight  on  D( [0,TJ ,$_q) . 


(i) :  The  log  characteristic  function  of  (x£  is: 


/  v'/.VAV  n'aV-V'' 


*•  w 


12 


0  ( f  •  •  •  ) 


[i  £  tk.k.»t»k)  + 


jo  jaxA  <eiaI?IP<8)1  -  1  -  iar^[P(s)])iin(dadr^)ds^j 


where  P(s)  J"  Vco,  tu]  <S)V 
k*l 


while  that  of  (Wfc  [$]cl)J,1  is 


C ( a«  ? •  • «  a— ) 


K  _ 

[i  V  tkakmltkl  “  ^  jo  rFC»>  >da  J  -  Hence 


lCn*al'  *  *  *aK*  ”  c<ai»  •  •  •  »aK^  I  *  I  *  /  tkak<m  ^k^  ”  m^k^  + 

k«l 


co 


lo  [La  L.  <iar?IF(s)l)V(dad^) 

P-3 


1  n 

— (Qn(F(s) /F(S) ) 
2 


Q(F(s),F<s)>) 


]]  ds|  < 


A 

i  y_  Vk‘""‘*iti  -  mt*ki|i ♦ 

Jc*l 

\ 

00 

lo  IexA  1  Y.  (iaI?lF(8)1)PlMn<dadI?,d8  * 

P-3 

fT  -|Qn(P(s),P(s))  -  Q(F(s),F(s)) |ds 

JO  A 


<  Y.  tkl#kl|B"[*lt1l  +  jo  jBxAl«^I»<s>l|3|in(dad^)as  + 


Jo  -|Qn<F<s),P<s))  -  Q<P(s) ,F(s) ) |ds 
2 

the  first  term  tends  to  zero  by  (5).  As  for  the  second  term,  use  (4) 
to  obtain 

(6)  lim  JrxA  larJlF(s)1 13^n(dadl7>  -  0  V  »  e  [0,T] . 

Now,  by  definition  of  F(s) 

K 

|a|^[P(s)]  |  <  Y_  lal  lakl  V  8  e  tO,T]. 

Define  a£  :■  a^  sign (a^  1^1^] ) .  Then 

i 

K  K  K 

Y  Ml»klll?l*kl|  ■  l»l  Y  »kI?l*kl  -  l»I?t  Y  “k^k'l 

ksl  k*l  k*l 


K 

(7)  |arj[P(s)]|  <  |ar^[  Y_  V  s  e  [o,ti. 

k-l 


I  a>k  «  i' 


so  an  application  of  (4)  gives 


But 


14 


lira  Jax/Jaljl  ^  a^k]  I3^n(dadl^>  -  0  and  thus 


!S£  /bxA111?1  Y.  •i»k>l3>,n(d,a,?)  <  “■ 

nera  k-1 


But  then  JbvA  |ap[F(s)  3 1  3pin(dad]^)ds - >  0 

n->oo 


by  (6),  (7)  and  the  DOT.  Further, 


Q  (FCs) ,F(s) ) - >  Q(F(s) ,F(s) ) 

n-»oo 


for  each  s  e  [0,T]  by  (3)  and  since  Q  and  Qn  satisfy  (1)  we  have 


|Qn(F(s),F(s))  —  Q(F(S)  ,F(s)  )  I  <  2C  H  F  ( S )  ||  ^  . 

r2 

Moreover,  H F C s )  ||  ^  ds  <  oo  so  the  DCT  gives 
jl  |Qn(F(s) ,F(s) )  -  Q(F(s) ,F(s) ) |ds - >  0 


concluding  the  proof  of  (i). 


n-*oo 


(ii)  By  (1)  and  Proposition  1.0  the  family  {p”  :  ne  ■>  is  tight 
on  D([0,T]  ,<Lq>. 


Let  A  :  |  |  b€  a  linear  and  continuous,  and  suppose  that  A  and 

{Tfc  :  t  £  0}  satisfies: 


There  exists  a  strongly  continuous  semigroup  [T^:  t>0] 


on  H  whose  generator  coincides  with  A  on  J  and  such  that: 

(a)  Tt$  cl  V  t>0 

<>»>  Tt|j  >  !  continuous  in  ($,  )  V  t  >  0 

(c)  t  is  "V  -continuous  for  every  $  e 

It  is  shown  in  [2]  (Theorem  III. 1.5)  that  the  SDE  on 
dxfc  -  A'Xfcdt  +  dMt?  XQ  -  Y, 

where  Y  is  a  -valued  random  variable  and  M  ■  < )  t> 0  is  a  $'“valued 
weak  L2-semimart ingale,  i.e.  a  -valued  process  such  that  Mt[$]  is  a 
semimartingale  with  E(Mt($l)2  <  oo  Vt>0  for  each  has  a  unique 

CADLAG  solution  on  . 

i 

Moreover,  the  following  result,  which  we  restate  for  the  convenience 
of  the  reader,  is  proved  ((2],  Theorem  III.2.1): 

Let  Mn  and  M  be  weak  $>' -valued  semimartingales  satisfying 

(*)  V  T>0  3  qT  e  «0  V  n  e  ■  :  Mn'T,  MT  6  D([0,T],$_q  ) 
and 

sup  E  sup  I|m£'T||2  <  oo. 

nea  o<t<T  -qT 

Let  Yn,  Y  be  -valued  random  variables  satisfying 


<**> 


-s.  ^  X.’'.  S.  ««-  ” 


q  r.  e  ■_  s  sub  max{E||Yn||ir  rE||Yfl2__  }  <  oo. 

J  1  °  new  rl  rl 

Let  T>0  and  suppose  that  Mn,T  *»>  MT  on  D([0,T],Ji  )  as  n-»oo 

hT 

and  that 


Yn  ■»>  y  as  n-»oo.  Then  there  exists  pT  e  «0  such  that  ^n,T  »»>  jr1 

on  D(  [0,T]  )  as  n->oo,  where 

PT 


:n,T  and  denote  the  unique  CADLAG  solutions  to 


-  A'jrjjdt  +  dMt;  To  “  Y°'  °-t-T 


respectively 


Tt  *  ^t'*w  ’  v“*t‘ 


.dt  +  dM.  ; 


Y,  0 < t <T . 


Moreover,  Theorem  1.2  remains  true  if  the  spaces  D([0,T],  §[_  J  and 


D(  [  0,  T] ,  <5 _ ]  are  replaced  by,  C(tO,T),  $_  1,  respectively  C([0,T] 


3L  1?  see  121 
C*T 


Let  and  J^°  be  ^'-valued  random  variables,  and  let  j£r  1  “  (J?'t>o 
denote  the  unique  solution  to  the  SDE  on  *§' 

*j?  -  +  si  -  s„ 

/ 

and  let  »  **?t*t>0  denote  tfle  unique  solution  to 


V  > ‘j’*’.-  v  v.sav>,.  A--.- . 


<JJ}t  -  A'l}tdt  +  <JWt;  -  T} 


-Remark  1 


If  the  operator  A  is  selfadjoint  and  dissipative  (regarded  as  a 
densely  defined  linear  operator  on  H)  with  <I-A)  being  a  Hilbert- 

Schmidt  operator  for  some  r1  >  0,  and  if  $  is  the  nuclear  space 

generator  by  (I-A)  (i.e.  $  -  6  Hs  ||(I-A)r^||g  <  co  V  r  e  R;  see 

[5])  then  the  solution  may  be  expanded  as  a  series 

co 

it  - 1  ^ 

(Converging  uniformly  on  10, T]  in  the  $_q  “  topology  (a.s.)  for 
every  T  >  0  and  q  >  rx  +  r2;  where  (fy  -\j);  0  <  <  \2  <  . . .  with 

\n  loo  as  n  oo,  is  the  eigensystem  for  A,  and  where  is  the 
one-dimensional  Ornstein-Uhlenbeck  process  given  by 

di?t  -  -s  ntat  +  • 

A  similar  expansion  is  possible  for  x£  in  this  case.  We  refer  to  [5] 
for  details. 

Remarft  g 

Regardless  of  the  structure  of  $  one  can  show  that  whenever  A 
satisfies  Al,  for  every  T  >  0  there  exists  PT  >  0  such  that  ( 


where  C([0,1],  i£“P^)  denotes  the  complete  metric  space  of  all 
continuous  functions  f:  [0,T]  <g  and  where 

-1*  (I?t,teio#T]  * 

1.3  THEOREM 


Suppose  that,  in  addition  to  (1), 


(8)  Qn t4»r4») - >  Q(4>,4>)  V  *  e  $ 

n-^oo 


(9)  lia»  J^I*<l?l*llV<d-dI}>  -  0  V  4>  e  $ 


(10)  3  rei!  sup  max{E||  ]^°j|^r,E||^n||_2>  <  co  and 

n 

■>  D°  on  $_r  as  n  4  oo. 

i 

i 

t 

(11)  mn[ 4»] - >  m[ 4>]  V  4  e 

n-»oo 


Then,  for  any  T  >  0  3 PT  e  ■ : 


r”'T  — > 
^  n->oo 


on  D((0,T],$ _ ) 

PT 


where  Jn'T  -  <£>tet0jT,  and 
r?T  “ 


PROOF: 


(1),  (8),  (9)  and  (10)  imply  that  Xn'T  —  ->  WT  on  D([0,T],<&  > 

n->oo  q 

V  T  >  0  by  theorem  1.1.  Moreover,  (1)  implies  condition  (*)  Theorem 
1.2  while  (10)  supplies  the  remaining  assumption  of  Theorem  1.2,  from 
which  the  conclusion  is  therefore  obtained. 

I 

Next,  we  shall  give  conditions  under  which  the  processes  Xn'T  will 

converge  weakly  on  D([0,T],$  )  to  a  process  XT  constructed  from  a 

Poisson  random  measure  N  on  RxAx[0,ao)  in  the  same  way  as  Xn  was 

constructed  from  Nn.  We  shall  then  invoke  theorem  1.2  to  give 

sufficient  conditions  for  the  weak  convergence  of  y  '  on 

D([0,T],<£  )  to  the  solution  to  the  SDE  driven  by  X. 

PT 

Let  m  e  <j>'  and  let  p  be  a  ^-finite  measure  on  (RxArB  (Mxfl(A) ) 
satisfying 

i 

(lia)  m[ 4>i 2  +  b(4>,4>)  <  c||<HI?  V  <j>  e  $ 

r2 

(llb)  IrxA  -  1  -  iaJ}l4>]|M(dadT})  <  oo 

where 

B(<M>  *-  a2JJ[^)]2M(dadrj);  V  <j>  e  $. 

Let  N  be  a  Poisson  random  measure  on  (RxAx[0,oo)  ,  fl  (Rxfl  ( A  )xfl  ( [  0,oo) ) ) 
with  intensity  measure  ^(dad]^)dt  (a  e  R,  6  Ar  t  >,  0). 


Define 

*t($)  *  j*xAxtO,T]  (N(dadI?da>  ~  M(dadrj)ds)  ?  t>  0;  <j>  6  j$, 

and  Xt(^)'  -  tmt$]  +  Yfc($). 

Since  the  ^  required  in  (11a)  is  the  same  as  that  of  (1),  theorem 
III. 1.12  (b)  of  [2]  implies  the  existence  of  a  $  valued 
semimartingale  X  ■  *xt*t>0  sat^s^y^n9  xtt4>l  “  Xfc($)  (a.s.)  for  every 
4  6  <$.  For  each  T  >  0  let  denote  the  measure  induced  on 
D([0,T),$__q)  by  XT  -  <xt>te[0,Tr  Th*n  bavft: 

1.4.  THEOREM 

Let  mn  and  jin  satisfy  (1).  Let  m,i*  satisfy  (llafb)  and  suppose  that 

\ 

U2)  '  /RxA  (e1*1?1*1  -  1  -  ia]^[4>]  )|in(dad]^) — -— > 

/■x A<eiaI^1*'  -  1  -  ia  1  )|i(dadj^)  V  <►  e  $ 

(13)  mnH>]  — »  ml*]  V  4  e  $. 


Then,  for  every  T  >  0,  p£  —>  as  n  4  oo. 


Fix  T  >  0.  Since  (1)  is  assumed  to  hold  {Pa  :  n  >  1}  is  tight  on 
D( 1 0,T] )  by  Proposition  1.0.  Hence  it  suffices  to  show  finite 
dimensional  convergence: 


Let  0  <  tx  <. ..<  tR  <  T  and  $k  6  <J>;  k  *  1,...,K. 

Then  the  characteristic  functions  for  (Xa'T[i(L  ], 

<xt  ^i  1  *  •  •  •  fXt  are'  respectively, 

X  K 

Cnul'***'aK)  *  exp|7mn[  ^  tkak^,k1  + 

ifci 


Jo  JrxA  <eiaJ?tF(s)1  -  1  -  ia]^[F(s)  ]  )pn(dad]^ )dsj 


.fXj,TnJrKJ)  and 
K 


C(alf...,ak)  -  exp |j.ml  WM  + 


Jo  JrxA  <eial^lF(S)1  "  1  -  iaI}[F(s)])n<dadr})ds 

K 

where  F(s)  :»  ak*[0,t  ]*s^k*  (13)  it  is  enough  to  show  that 

k=l  '  k 

n^oo  CXP[io  l»xA(*ia^[F(a>1  “  1  “  i*IJ[F(s)]  )Mn(dadr))dsj  - 
*Xp[]o  JrxA  (eiaI?lF(s)1  -  1  -  ia D [F  ( s)  ] p  (dadD  ) ds]  . 


Now,  F ( s )  is  piecewise  constant,  i.e.  there  are  0  »  s 
and  e  $  such  that 


F  ( s )  = 


f  $ j  if  S  e  CSj_1,Sj)  j  »  1, . . .  ,M  -  1 
WM  if  S  6  [sM_1  ,T3 


Hence 


eia]}[F<s)l  x  _ 


ia  J^[F(s)  ] 


M— 1 


)  le  3  -  1  -  iamf.]>l[s  ,(s)  + 

3^=1  '  3 


(e 


1  -  iaI}lV,ltsM_1(T]<a' 


SO  ; 


lo  JrxA  (eial?tF(s)1  -  1  -  ia^[F(s)]Mn(dadr^)ds 


M— 1 


*xA 


(e 


ia  14»  j  1 


ia]^[^j]  )^n(dadl^) 


lr.  _  v(s)  +  |RxA  (e 


Uj-l,aj) 


ial^l^l 


-  1 


-  l.lJ[*h])Mn(d.dJJ)l[  ,  ,T)<s>"|da 


2*A 


V 

L  'RxA 

3-1 


(e 


iaj|  L  9j  i 


3  _ 


1  —  ial^[^j  ]  )(jin(dad]^)  (Sj  —  Sj^) 


(by  (12) 


M 


— >  V"  [  (6ia^[+j1 

_  .i  /  JrxA  ( 


n-»oo 


3=1 


-  1  -  ia]^t<^j])ji(dad]^)  (Sj  — 


(Recall  that  |RX//^  ...p(dadl^)  is  finite  by  (lib)) 
“lo  jaxA  (eiaI?[F(s)1  -  1  -  iaJ}[F(s)]  )n(dadl})ds. 


concluding  the  proof. 


} 


5°  be  a  <£'— valued  random  variable  and  let  £ =  denote 

the  unique  solution  to  the  <£'—  valued  SDE 


dlt  -  A'5tdt  +  dxt'  Jo85!' 


1.5.  THEOREM 

Let  m11  and  pn  satisfy  (1),  let  m,^  satisfy  (11a, b)  and  suppose  that 

(12)  and  (13)  hold.  Suppose  further  that 

(14)  ]  r  6  B  :  sup  n»ax  {E  |  lr>E  |  5°|lir>  <  ao 

and  that  J  »=»>  J°  on  . 

n-»oo  -r 


Then,  for  any  T  >  0,  3  pT  e  H  : 


■“>  on  D( [ 0,T] ,(E  ),  where 

n->ao  PT 


(5t)tetO,T]. 


PROOF: 


Let  T  >  0.  Recall  that  q  >  r,  is  such  that  the  canonical  injection 
*  r2 

■Cg  is  Hilbert— Schmidt  from  <£g  (£r  •  Let  <4> j  :  j  6  ■>  be  a  CONS 

2 

in  <$q  consisting  of  elements  of  <$.  Then  note  that 


ao 


E  sup  ||xJ'T||2  «  E(sup  V  (xJ,T[«j>.])2)  < 

0<t<T  c  q  0<t<T  i—  3 


CD 

E(7  sup  <x"'T[4>.1)2> 
0<t<T  c  3 


ao 

Y  Esup  (x“[<|>.])2  < 
0<t<T  11  3 


oo 

y  2(T2mn[<J>j]2  +  4TQnt|j(^|)  <  (by  <1>) 

j-1 

i 

ao 

Y  2C(T2v4T)  8 ^  =>  ZCer^DlK  2|2g  V  nB  B, 

pi  2 

(where  II  *  It HS  denotes  the  Hilbert-Schmidt  norm)  i.e.  (*)  of  Theorem 

1.2  is  satisfied.  Moreover,  Xn'T,XT  €  D([0,T],<5>  (P-a.s.)  by 

'  q 

assumption  and  Xn  and  X  are  <£'  —  valued  (weak)  L2— semimartingales.  By 

Theorem  1.4,  (1),  (11a, b),  (12)  and  (13)  imply  that  Xn'T  —->  XT  on 

n->oo 

D( 10,T] ,!L_q)  •  Since  also  (14)  is  supposed  to  hold,  the  assumptions 
of  Theorem  1.2  are  satisfied  and  the  conclusion  therefore  follows 


from  this  theorem. 


-Next  we  shall  give  conditions  for  the  weak  convergence  of  a  sequence 
Wn  of  <$'— valued  Wiener  processes  to  another  <$'— valued  Wiener  process 
W r  and  then  employ  these  together  with  Theorem  1.2  to  giving  the 
corresponding  weak  convergence  result  for  the  solutions  to  the  SDE's 
driven  by  Wn  and  W,  respectively. 

Let,  for  new,  mn  6$'  and  let  Bn  :  <|x$  i  be  bilinear  symmetric 
functionals  satisfying  (1).  Let  Wn  ■  (wtJt>0  denote  the  —valued 
Wiener  process  with  parameters  m11  and  Bn.  (1)  implies  that  w“  6  $ 

V  t  >  0,  for  some  q  which  does  not  depend  on  n  6  H. 

1.6.  THEOREM 

Suppose  that,  in  addition  to  satisfying  (1),  Bn  and  mn  satisfy 

(15)  Bn(^>,<|>) - >  Q(<|>,<fr>  V  M  1 

j  n-*ao 

(16)  mn[<fr] - >  m[<j>]  V  ^  e 

n->ao 

Then,  for  each  T  >  0,  we  have  Wn'T  ■■■>  wT  on  C(  [0,T]  ,<E  ) , 

n-»ao 

where  Wn'T  »  *wt*t£LQ,T]  and  wt  ia  the  —valued  Wiener  process 
introduced  on  page  6. 


PROOF: 


By  Mitoma,  [6]  (Theorem  5.3  part  1),  it  is  sufficient  to  prove  that 


V  t  >  0  {W11'1  s  n  6  1}  is  tight  on  C(tO,T],<5 _ )  and 

4 


V  o  <  tx  <...<tN  <  t  V  e  It 


<w?;T  <wV^ni-i- 


The  tightness  part  is  proved  in  the  same  way  as  in  Proposition  1.0 


Now,  a  calculation  shows  that 


N 


E  exp( i  Y"  a.,w“'T  [f.])  - 

h  j  j 


N 


6XP[i  tjaj®ni^j1 

1-1 


.  N  N 

A  V  V 


“II 


j-1  k-1 


(by  (15)  and  (16)) 


n-»ao 


N 


_  !  JL  JL 

exp[l  --  X^jVjVW] 


j-1 


j-1  k-1 


N 


E  exp(i  \  a.W^  [f,]  ). 


k^ 


} 


Letting  rf n  ■  ( t  ^  t > 0  denote  the  unique  solution  to  the  SDE  on  $' 


d^“  -  A'^dt  +  dwj  ; 


be  the  J'— valued  process  introduced  on  page  16  we 


and  .  <r}t>t>0 
have 

1.7.  THEOREM 

Let,  in  addition  to  (1).  Bn  and *ran  satisfy  (15)  and  (16)  of 
theorem  1.6,  and  suppose  that  ]^n  and  ]^°  satisfy 


(17)  3  r  6  *  *  auP  max{E||  J}n||lr'  Elf}°ll-r*  <  00  and 

n 


IJn  ~ on  $-r' 


n-»oo 


Then,  V  T  >  0  3  pT  6  ■:  I)n'T  ~>  r>T  on  C([0,T],»  ), 

*  n-»a> 


where  I}n'T  v 


<r}t}teto,T] 


1 

PROOF : 


Bv  (1),  (15)  and  (16)  and  theorem  1.6  Wn,T  — >  WT  on 

•r,  n->ao 

D(tO,T],<$_  )  V  t  >  0  where  q  ■  min{p  :  «£_  is  Hilbert-Schmidt ) . 
Moreover  (1)  implies  (*)  of  Theorem  1.2  and  (17)  supplies  the 
remaining  condition  of  Theorem  1.2  (recall  remark  2). 


As  indicated  at  the  beginning  of  this  section,  Kallianpur  and  Wolpert 
(15)]  used  Poisson  random  measures  fie/ined  via  intensity  measures  on 
(■x  X ,  8  (M)xfl)  where  is  a  suitably  chosen  measurable  space, 

rather  than  by  mean/covariance  measures  defined  on  (IxA» fl(R)xB(A> )  • 


A  €  a($' )  as  we  have  done  it  here. 

It  is  therefore  natural  to  address  the  question  of  when  Kallianpur 
and  Wolpert's  framework  is  contained  in  the  one  we  have  presented 
here. 

It  is  proved  in  [2]  (Chapter  IV)  that  if  H  -  L2(X,T)  for  some  - 
finite  measure  space  (X,T)  where  X  is  a  -compact  topological 
Hausdorff  space  then  it  is  possible  to  represent  the  Kallianpur  & 
Wolpert  models  in  our  framework  if  elements  of  ^  are  continuous 
functions  on  X  and  if,  the  evaluation  functional  &x  :  $  $(x)  is 

continuous  on  $  for  every  x  e  X  • 

II.  THE  WAN  &  TPCKWELL  MODEL 


Next,  we  shall  apply  our  results  to  giving  a  rigorous  formulation  and 
investigation  of  a  model  recently  proposed  by  Wan  6  Tuckwell  [91s 

/ 

In  order  to  study  the  behaviour  of  the  difference  V(t,x)  at  time  t 
between  the  so-called  resting  potential  and  the  actual  potential  at 
point  x  on  the  surface  of  an  infinitely  thin  cylinder  shaped  neuron 
which  receives  synaptic  stimuli  of  the  finite  spatial  extent  6^  at 
each  of  N  sites  x^,  Wan  &  Tuckwell  investigated  the  model  formally 
given  by 


iV 

/"  - —  X 

J2V 

_N_ 

^  h(x;xi,ei) (W\  + 
i-1 

dWA 

A — > 

1  dt 

1 

(19) 

■  V  T 

it 

,*2  + 

i 

S  V(0,x)  -  0 

vct,o> 

-  0  -  V( t, b) ;  V 

t  >  0, 

where 
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-  i<xl-«i.xl*€1>u' 

(x^,e^  >  0  fixed  for  i  -  1,...,N) 

and  where  w*;  i  »  1,...,N  are  independent  standard  Wiener  processes. 

and ^4^  represent  input  current  parameters  and  the  neuron  is 
thought  of  as  the  interval  (0»bl;  for  some  b  >  0. 


To  see  how  this  model  can  be  given  a  rigorous  representation  as  a 
valued  SDE,  let  H  -  L2(I0,b])  with  inner  product  denoted  by 

Let  L  denote  the  operator  I  -  A  <A  *  Laplace  operator  in 
one  dimension)  with  Neumann  boundary  conditions  at  0  and  b.  Then  L  is 
a  densely  defined  positive  definite  selfadjoint  closed  linear 
operator  on  H  and  admits  a  CONS  :  j  ■  0,1,2,...}  in  H  consisting 
of  eigenvectors  of  L; 


j  “  0,1,2,...,  where 


1  + 


and 


<frj(X)  - 


r  b“1//2  if  j  »  0 


jlfx 

Cos  ( - )  if  j  >  1. 

b 


Further,  A  :«  — L  is  the  generator  of  a  selfadjoint  contraction 
semigroup  {Tfc  :  t  £  0}  on  H  whose  resolvent  R<\)  ■  (\I  —  A)  *  is 
Hilbert— Schmidt  on  H. 


Letting 


5^:-  <<fr  e  H  :  I  (I  -  A)r^flH  <  oo  V  r  e  R} 
and  defining  norms  fl.|r;  r  6  R  on  $  by 

II  $  II  r  |  (I  -  A)r<|>||H;  4>  6  $ 

we  put  <$r  equal  to  the  | .  |  r~completion  of  J. 

Then  $  ■  and  if  V denotes  the  Frechet  topology  on  $  generated 

by  { H  .  Q  :  r  e  R>  (i.e.  the  projective  limit  topology  on  $}),  then 
($,I  )£•>  H  (where  denotes  the  strong  dual  of  ($,0 )  is  a 

rigged  Hilbert  space.  Since  A  ■  — L,  and  L  is  a  densely  defined 
positive  selfadjoint  closed  linear  operator  on  H  it  is  easily  seen 

that  A  and  {Tfc  s  t  >  0}  satisfy  A1  of  section  I. 

\ 

i 

Moreover,  (^  :  j  6  1}  I  Dcm(L)  and  per  construction  of  $ 

every  element  of  <$  is  an  infinitely  differentiable  function.  Let 
N  €  H  fixed,  and  for  each  i  «  1,...,N  let  Ji  el', 
i  -  1 , . . . ,  N  be  f*-finite  measures  on  R  satisfying 

ja  a2vi(da)  <oo  Vi, 

and  let  \i  be  the  measure  on  RxA#  where  A  ■  (J^  s  i  ■  1»...#N}» 
given  by 
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7  ;  where  is  the  point  mass  at 

U  5 


Define 

Q<<M>  -  JExA  a2r^t^ir|n|r]M(<aad^);  4»r4i  e  i 

N 

-  iu  aS(da)  iii^?i1^ 

then  Q  is  a  continuous,  bilinear  symmetric  functional  on  $,  so  for 
m  e  given,  let  W  ■  Wfc  be  the  <|j' -valued  (actually  $  valued  for 
some  q  e  flQ)  Wiener  process  with  parameters  m  and  Q. 

Consider  the  SDE  on  <$' : 


(20)  dl}t  -  A'I}tdt  +  dWt,  I}0  -  0 

*i 

I  y 

Now,  w  is  a  weak  (£'— valued  continuous  L  —semimartingale,  and  since  A 
and  (Tt  :  t  ^  0}  satisfy  A1  there  is  a  unique  continuous  valued 
solution  (from  [2],  Theorem  III. 1.5  and  Remark  6)  given  by 

-Jo  Ws[Tt_sA^]ds  +  Wt(*]  VMI  (with  probability  one). 

N 

Choosi  119  V  ^  *  If  •  •  •  f  H  &nd 

N  , 

®  —  ®  *  —  'y  ^ |  ^ m  Jr  *  ■Wj  (da ) , 

(20)  is  the  representation  of  (19)  as  an  SDE  on  .  To  see  that  this 
is  indeed  the  case,  expand 


4> 


/  (converging  in  ($,D) 

3-0 


j recall  that  ^  e  |  V  j  6  ■) 

Then  (writing  for  and  W®  for  Wfc) 
ao 

IJ®t+l  -  Y  <£  »fl*t-sA*jIds  +  “?'♦  j>>  <*-*3>H 

J=J 

2 

(converging  in  L  (&,F,P)). 

Moreover,  the  series 
oo 


Y  Oo  tTt-,A+j ,d* ♦  ^♦ji]+j<«» 

is  convergent  in  L2(A,P,P)  for  every  x  e  [0,b]  to 

satisfying 

; 

^t1*1  *  jo  Ve(t,x)^»(x)dx  (P-a.s.) 


for  every  $  6  <£.  Let  ns  see  that  EVg(t,x)  and  VarVg 
agree  with  the  formulae  found  in  [9]  by  a  heuristic 


EVg ( t , x ) 


E 


wsITt-sA*j,dS  +  wtl*J,,*j(x) 


«  E 


(t-s) 


limit  Ve(tfx) 


t,x)  actually 
argument : 


W_H>,lds  +  W  [«!>.] )4> .  (x> 


dWsf$j (x) 


■  T  II 


-X4(t-s) 


00 


1=0 


(x) 


ao 

Y  ■l*jlVj1(l 
J=0 


-Xit 

e  -1  )^j(x) 


f  ^  f  ti^Wa  .  .’V, 


£i  1  £o  xj 


which  is  fornula  (8)  page  279  in  Wan  6  Tuckwell  19] .  Here,  as  in  19] 


|j(xi;€i)  -  <h(.;xi,ei)^j>H  »  (  1  **  4>j 


(x)dx. 


Next, 


ao  ao 


VarVg  ( t 


f  f  -Xj (t-s)  . 
'X)  “  JO  e  3 


-Xv<t-s)  . 


dW9(^ltJ^j  (x)$k(x) 


TT  ft  “*^i+XjtHt-s> 

Jo  *  3  Q(^j,^k)ds^j(x)^k{x) 


j^)  )c£& 


v  ^2,Q(+i^Jc)1  "(Xi+Xk)t 

/  /  - 3-JL-^  (x)a  (x>  (i  -  e  3  *  ) 


f.2f  v-*i(x,V*>V*i’ei>V*i’ei> 

l  l.  r~n - 

i-1  j-0  k»0  Kj  *  Kk 


oo  ao 


(1  -  e  3  * 


) 


which  is  formula  (10)  in  [9] 


Wan  &  Tuckwell  proceed  to  compute  the  limit  as  6^  -»  0  V  i  *  1,...,N 
in  such  a  way  that  ->  and  €i/(^ i  -»  >  0  of  EVg(t,x)  and 

VarVg(t,x),  and  they  find  that  these  limits  correspond  to  having 
point  stimuli  (i.e.  h(x,xi,ei)  replaced  by  £x  (x))  at  each  of  x^ 


1 1  •  •  •  t  N. 


This  result  may  be  obtained  from  Theorem  1.7  in  the  following  manner: 

For  each  i  »  1,...,N  take  »  b^ej1  »  where  a  is  finite 
measure  on  R  with  compact  support. 

Noting  that  every  $  6  is  a  continuous  function  on  [0»b]  (recall 
that  <$CDom(L)  and  that  L  is  a  differential  operator)  we  let  6^  -»•  0 
in  such  a  way  that  ->  a^.  Then 

N 

V* tfK*.  <h( .  ;x.  ,e.  )  ,4»_ 


lim  mg[$]  -  lim 


and 


1 24i^i 

i»l  3 


[<j>] 


lim  Q6  ( 4>  r  <t> ) 

ei->o 


N 

Um  Ya\ 


(Xi+€i<j>(x)dx)2 

Jxi-Si 


im  Y  b?e72(  I*1  6i^(x)dx)2 

Jxi-Si  * 


lim 

e 


(da) 


^T«>i*‘V2 


r.» 

J_  “i 


(da) 


Mi (da) 


Also, 

\ 

f  mec4»]  {  2  +  Q6(<f>,<!>)  < 


H 

[</"  |h«.i*i.«1)|Hl^l|)2  + 
f=i 


N 


T/Ub< 

i-l 


'ei> 


«  ( 


N 

y  2«i i  <^i i 


+ 


N 


2 

H 


<  CONSTANT  H  4>  II  h  V  «i» 


since  6^  -»  ai  and  eL  -»  0;  where  CONSTANT  is  independent  of  so 
condition  (1)  of  section  1  is  satisfied.  Since  the  initial  condition 
is  zero,  theorem  1.7  yields 


Proposition 

As  ei  -»  0  and  e^ac^  -»  a^  we  have 


®lT  «P 

T}  1  — =>  j}1  on 


C( [0,T] ,  §>  )  V  T  >  0 

“qT 


for  some  qT  >  0,  where  1^  *  i-8  the  solution  to  (20) 

corresponding  to 


_N_ 

Q(<|>,4>)  -  )  4  b?t<|)l 

frl 


^(da) 


and  ; 

N 

m[<t>]  =  ^  2a.  £x  I*]  . 
i*L  1 

Now,  Take  R  a2  Mi(da)  *  1  i*l,...,N.  Then 


N 

V 
L . 

i-1 


2a  i/__  ~ - 

0 


(1 


j* 


♦  6  $ 


and 


N 


Var 


rjtt+i  -  y  4b iy  y  ^'Vg^'Vgf^v+fcQ,) 

j  ,  r-»  v  .  v 


j  T  *k 


and 

V-  2  V5-  ^-4  (xi ^v(xi * 

(23)  VarVg( t,x)  -  /_  4bi  2_  — ~ — “ — — <t>j  (x)^k(x) . 

i»i  j-0  k=0  Xj  +  Xk 

(22)  and  (23)  are  the  expressions  found  by  Wan  &  Tuckwell  for  point 
stimuli  at  x^i  »  1,...,N. 

In  practice,  equation  (20)  is  likely  to  arise  as  a  limit  of  equations 
where  the  noise  is  not  u  Wiener  process,  but  rather  a  process 
generated  by  a  Poisson  random  measure  in  the  manner  considered  in 
section  I.  As  an  illustration,  take  nn  to  be  measures  on  RxA?  where 


where 


for  each  n  6  H  and  i  *  is  a  C-—  finite  measure  on  a  such 

that 

sup  f_a2**?(da>  <  C  <  ao  Vi*  1,...,N. 
n€H  K  1 

Let  m11  6  converge  weakly  to  nr€.  Then  there  is  r  6  HQ  such  that 

<  K|*nJ  V  n  S  *. 

Since 

fJiH>H2  <  <2ei)2||<|»|| ^  <  (2ei)2j^||2  we  get 

— 

|  mn  [  <f>  ]  |  2  +  Qn(<t>,<j>)  -  |mn[$]  |  2  +  )  jRa2-V?(da)  ( ^[^D2 
<  CONSTANT  i|<|>||2  y  n  e  H;  i.e.  (1)  holds  with  r2  -  r. 

i 

Let  x”;  n  >L  1  denote  the  (£' -valued  processes  constructed  from  mn  and 
Un  on  p  .S'/i  . 

Letting  denote  the  solution  to 

\ 

dfj  -  -L'jj?  +  dXt 

SS-0 


Theorem  1.3  gives  the  existence  of  pT  such  that 
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V*'T  — >  ne'T  on  D([0,T],$_  ) 

5  n->ao  (  pT 

provided  that 

(24)  lim  JR|a|3-4j(da)  -  0  Vi-  1,...,N 
n-*a> 

and 

(25)  lim  La2V?(da)  ~A  ?  V  i  - 

n->ao  **  1  '  1 

a 

i.e.  the  previously  considered  process  can  be  thought  of  as  the 
limit  of  solutions  to  SDE's  with  Poisson  generated  noise. 

Physically,  this  type  of  weak  convergence  models  a  situation  in  which 
the  individual  current  stimuli  of  the  neuron  arrive  very  densely  in 
each  small  time  interval  so  as  to  create  a  total  contribution  to  the 
electrical  potential  which  behaves  like  the  continuous  Wiener  process. 

On  the  other  hand,  if  (24)  and  (25)  are  replaced  by 

lim  la  (eiay  ”  1  ~  iay)-vj(da)  -  jR  (eiay  -  1  -  iayW?(da) 
n-»ao 

for  all  yea,  then  theorem  1.5  gives 

— «>  *e'T  on  D(C0,T],J  ) 

*  n-»co  *  Pt 

where  is  the  process  with  mean  functional  m®  constructed  from  the 
Poisson  random  measure  with  intensity 


This  latter  convergence  can  be  thought  of  as  modelling  a  situation  in 
which  the  individual  stimuli  received  by  the  neuron  do  not  tend  to 
arrive  very  densely  packed  in  each  small  time  interval,  but  rather 
tend  to  arrive  clustered  at  random  points  of  time. 
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